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Abstract. A fc-symplectic framework for classical field theories subject to 
OO , nonholonomic constraints is presented. If the constrained problem is regular 

C D . one can construct a projection operator such that the solutions of the con- 

C D ' strained problem are obtained by projecting the solutions of the free problem. 

' Symmetries for the nonholonomic system are introduced and we show that 

; ^ ^ for every such symmetry, there exist a nonholonomic momentum equation. 

, The proposed formalism permits to introduce in a simple way many tools of 

nonholonomic mechanics to nonholonomic field theories. 
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1. Introduction 

During the past decades, much effort has been devoted to the differential geo- 
metric treatment of mechanical systems subject to nonholonomic constraints. To 
a large extent the growing interest in this field has been stimulated by its close 
connection to problems in control theory (see, for instance, iZjilO ). In the litera- 
ture, one can distinguish mainly two different approaches in the study of systems 

Key words and phrases. Nonholonomic constraints, classical field theories, fc-symplectic for- 
malism, nonholonomic momentum map. 
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subjected to a nonholonomic constraints. The first one, commonly called nonholo- 
nomic mechanics, is based on the d'Alembert's principle. This principle specifies 
from the constraints a subbundle of the tangent bundle, representing the admissible 
infinitesimal virtual displacements. The second one is a constrained variational ap- 
proach called vakonomic mechanics [2]. As is well know, the dynamical equations 
generated by both approaches are in general not equivalent [11] . 

In this paper we will study an extension of nonholonomic mechanics to classical 
field theories with external constraints. Nonholonomically constrained field theories 
have already been studied in the literature. The mathematical framework for a 
nonholonomic field theory that has been proposed in fSJ involves, among others, a 
generalization of d'Alembert's principle and of the so-called Chetaev rule that is 
commonly used in nonholonomic mechanics to characterize the bundle of constraint 
forms representing the admissible reaction forces. The constrained field equations 
for classical field theories, are then derived in a finite-dimensional multisymplectic 
setting. In [12] the authors continue and extend the work described in [3] . 

The multisymplectic formalism, was developed by Tulczyjews school in Warsaw 
(see, for instance, [2V), and independently by Garcia and Perez- Rendon [T51 [H] 
and Goldschmidt and Sternberg |16j . This approach was revised, between others, 
by Martin [30l[3T] and Gotay et al [17] and, more recently, by Cantrijn et al [9] (see 
also [3S] and references therein). 

An alternative way to derive certain types of the field equations is to use the 
fc-symplectic formalism. The fc-symplectic formalism is the generalization to field 
theories of the standard symplectic formalism in Mechanics, which is the geomet- 
ric framework for describing autonomous dynamical systems. In this sense, the 
fc-symplectic formalism is used to give a geometric description of certain kinds of 
field theories: in a local description, those theories whose Lagrangian does not de- 
pend on the base coordinates, denoted by {t^, . . . , t'^) (in many of the cases defining 
the space-time coordinates); that is, the fc-symplectic formalism is only valid for 
Lagrangians L{q^, w^) and Hamiltonians H{q^,pf) that depend on the field coordi- 
nates and on the partial derivatives of the field w^, or the corresponding moment 
pf. A more general approach has been given in 29J using the fc-cosymplectic for- 
malism. 

Giinther's paper [19j gave a geometric Hamiltonian formalism for field theories. 
The crucial device is the introduction of a vector-valued generalization of a sym- 
plectic form, called a polysymplectic form. One of the advantages of this formalism 
is that one only needs the tangent and cotangent bundle of a manifold to develop it. 
In [51] Giinther's formalism has been revised and clarified. It has been shown that 
the polysymplectic structures used by Giinther to develop his formalism could be 
replaced by the fc-symplectic structures defined by Awane [S] Sj . So this formalism 
is also called fc-symplectic formalism (see also [^ \T7[ 

Let us remark here that the polysymplectic formalism developed by Sardanashvily 
[ini [Ml SO] , based on a vector- valued form defined on some associated fiber bundle, 
is a different description of classical field theories of first order than the polysym- 
plectic (or fc-symplectic) formalism proposed by Giinther (see also [20] for more 
details) . We must also remark that the soldering form on the linear frames bundle 
is a polysymplectic form, and its study and applications to field theory, constitute 
the n-symplectic geometry developed by L. K. Norris in [32 ] [35 ] [36 ] [37 ] [38] . 

The purpose of this paper is to give a fc-symplectic setting for first-order classical 
field theories subject to nonholonomic constraints. In the fc-symplectic setting 
we will construct, under an appropriate additional condition, a kind of projection 
operator that maps solutions of the free problem into solutions of the constrained 
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problem. Nonholonomic symmetries are introduced and we show that for every 
such symmetry, there exist a nonholonomic momentum equation which reduces to 
a conservation law when the constraints are absent. 

We analyze some particular cases, for instance, the case of a constraint subman- 
ifold M obtained as A:-copies of a distribution D in Q has special interest. In this 
particular case, we construct a distribution H on T^Q = TQ® . (BTQ (i.e. the 
Whitney sum of k copies of TQ) along M such that for each Wq ^ M, H^^ is 
a fc-symplectic subspace of the fc-symplectic vector space {Tw^{TlQ),uo\{wq), . . . , 
u)'l{wq); V{'Wq)) where {uj}^, . . . ,uj'1;V) is the /c-symplectic structure obtained from 
L. Thus if we restrict the fc-symplectic structure of T^Q to H, the equations of 
the constrained problem take the usual form for a free problem at each fibre of H . 
This procedure extends that by Bates and Sniatycki [5] for the linear case. 

The scheme of the paper is as follows. In Section [5] we recall some basic ele- 
ments from the fc-symplectic approach to (unconstrained) Lagrangian classical field 
theories. In SectionOwe discuss the construction of a nonholonomic model for first- 
order Lagrangian Classical field theories with external constraints and we obtain 
the corresponding nonholonomic field equations. Next, in Section |4] we construct, 
under an appropriate additional condition, a projection operator which maps so- 
lutions of the free problem into solutions of the constrained problem. In Section 
[5] we derive the nonholonomic momentum equation. In Section [H] we analyze some 
particular cases and in Section [7] we briefiy analyze the Hamiltonian case. Finally 
m Section El we conclude with some general comments. 

All manifolds are real, paracompact, connected and C°°. All maps are C°°. Sum 
over crossed repeated indices is understood. 

2. fc-SYMPLECTIC Lagrangian field theory 
2.1. Geometric elements. 



2.1.1. The tangent bundle of -velocities of a manifold. Let tq : TQ Q he the 
tangent bundle of Q. Let us denote by T^Q the Whitney sum TQ® . ®TQ of fc 
copies of TQ, with projection r : T^Q Q, t{vi^, . . . , Vkq) — q, where va^ £ TqQ, 
l<A<k. 

T^Q can be identified with the manifold Jq (M'^, Q) of the k^ -velocities of Q, that 
is, 1-jets of maps cr : M*^ ^ Q with source at G M*^, say 

Jl(R^Q) = TQ®.^.®TQ 

J0,qC^ = (t'lq,---,'yfcq) 



d 

where q = cr(0), and fAq = cr*(0)(^p- 



). TlQ \s called the tangent bundle of 



k^ -velocities of Q (see [33]). 

If (q*) are local coordinates on [/ C Q then the induced local coordinates (g*, w*), 
1 < i < 71, on TU — Tq^{U) are given by 

and the induced local coordinates (g*, v\), 1 < i < n, 1 < A < k, on T^U = t^^{U) 
are given by 

g'(wiq, . . . , Vkc^) = q'{q), v\{Vlc^, Wfeq) = VAq{q') ■ 
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A. Vertical lifts of vector fields from Q to T^Q. 

Definition 2.1. For a vector X^^ G TqQ, and for A — 1, . . . , k, we define its vertical 
A-lift (Xq)^-^ as the local vector field on the fiber T^^(q) C T^Q given by 



{Xq)^'^{'Wq) ^ ^(«lqj • ■ ■,VA-lq,VAq + sXq,VA+l^ 

for all points Wq = (uiq, . . . , Vk^) G t^Hq) ^ T^Q- 

d 



,Vkq) 



s=0 



In local coordinates, if X„ = a* . 



then 



d 



(2.1) 



If X is a vector field on Q then we define its vertical A-lift to T^Q, 1 < A < fc, 
as the vector field X^^ given by 

d 



(X*or)K) 



then 



where X = X^ 



X"^^ = {X' o r) 



A 

_d_ 



B. Complete lift of vector fields from Q to T^Q- 

Let $ : Q Q be a differentiable map then the induced map T^^ : T^Q T^Q 
defined by T^^{j}^a) = jo(<f> o a) is called the canonical prolongation of $ and, it 
is also given by 

where -wiq, . . . , Ufeq G TqQ, q G Q. 

If Z is a vector field on Q, with local 1-parametric group of transformations hs : 
Q ^ Q then the local 1-parametric group of transformations T^{hs) : T^Q — > 
which is the flow of the vector field Z'~' on T^Q, called the complete lift of Z. Its 
local expression is 



d_ 



dZ'' d 
' dq^ dv\ 



(2.2) 



where Z — Z^-—. 

dq' 



C. Canonical fc-tangent structure. 



The canonical k-tangent structure on T^Q is the set {S^ , . . . , S'^) of tensor fields 
of type (1,1) defined by 

5^(«;q)(Z^J = (r,(«;q)(Z^J)^^(w;q), for all G T^^ (T^^Q), w^eT^Q, 

for each A — 1, . . . ,k. 

From (|2.ip we have in local coordinates 
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A 

The tensors S can be regarded as the (0, . . . , 0, 1 , 0, . . . , 0)-hft of the identity tensor 
on Q ioTlQ defined in Morimoto [55] , 

In the case k = 1, is the well-known canonical tangent structure (also called 
vertical endomorphism) of the tangent bundle (see ITS l [22]). 

D. Canonical vector fields. 

Let us denote by A the canonical vector field (Liouville vector field) of the 
vector bundle t : T^Q ^ Q . This vector field A is the infinitesimal generator of 
the following flow 

Tp-.M.XTlQ >TIQ , V(s,'flq, • ■ • j-ffcq) = (e^Wlq, . . . jB^Wfeq) , 

and in local coordinates it has the form 

n k t~. 



i=l B=l 



B 



A is a sum of vector fields Ai + . . . + A^, where each A^ is the infinitesimal 
lerator of the following flow ip"^ -.RxT^Q — > T^Q: 



1p {s,Vlq, . . .,Vkq) = (Wlq, • ■ . , V VAq, V A+lq, ■ ■ ■ , ffeq) 

and in local coordinates each Aa has the form 

AA^y"v\ — ^ l<A<k. (2.4) 

2.1.2. Second- order partial differential equations in T^Q- 
fc-vector fields and integral sections. 

Let M be an arbitrary smooth manifold. 

Definition 2.2. A section X : M — > T^M of the projection r will be called a 
/c- vector field on M. 

Since T^M is the Whitney sum TM® .'f . ®TM of k copies of TM, we deduce 
that to give a fc-vector field X is equivalent to give a family of k vector fields 
Xi,. . . ,Xk on M defined by projection onto each factor. For this reason we will 
denote a fc-vector field by {Xi, . . . , Xk). 

Definition 2.3. An integral section of the k-vector field X — {Xi, . . . , Xk), passing 
through a point q G M , is a map ■0 : C/o C M*-^ ^ M, defined on some neighborhood 
Uo ofOeR'", such that 



<')(J.|,) 



V'(O) = q, ^,(t) — ]^XAm)) , fortef/o, l<A<fc 



or, what is equivalent, satisfies that X o ip = ip'^^\ where ^^^^ is the first prolon- 
gation of ip to T^M defined by 



A k-vector field X — {Xi, . . . , Xk) on M is integrable if there is an integral section 
passing through every point of M . 
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In local coordinates, we have 

ij^'\t\...,t'^)=(^i;\t\...,t'^),^{t\...,t'')^ ,l<A<k,l<t<n, (2.5) 

and ijj is an integral section of (^i, . . . ,Xk) if and only if the following equations 
holds: 

J^^{XAyo4, l<A<k, l<i<n. 

In the fc-symplectic formalism, the solutions of the field equations are described as 
integral sections of some fc-vector fields. Observe that, in case A; = 1, the definition 
of integral section coincides with the usual definition of integral curve of a vector 
field. 

Second-order partial differential equations in T^Q. 

The aim of this subsection is to characterize the integrable fc-vector fields on 
T^Q whose integral sections are first prolongations (/)*^^^ of maps : R'' Q. 

Definition 2.4. A k-vector field ~ (^i, . . . , £^k) on T^Q, is called a second order 
partial differential equation (sopde) if it is a section of the vector bundle T^t : 
Tl{TlQ)^TlQ;that is, 

Tkr0{^l,...,^k) = IdrplQ , 

or equivalently 

r*(wq)(CA(«^q)) = for all A^\,...,k, 

where = (wi^, . . . ,Vk^) G T^Q. 

In the case k — 1, this is just the definition of a second order differential equation 
fSODEj. 

From a direct computation in local coordinates we obtain that the local expres- 
sion of a SOPDE ^ = (Cl, • • ■ I'C/c) is 

U{q\v\)^v\^^+{UrB^. l<A<fc, (2.6) 
oq avg 

where [S^aYb ^'"'^ functions on T^Q. 

If : ^ T'fcQ is an integral section of ^ = (^i,...,^^), locally given by 
'0(t) = (V'*(t), -(/^^(t)), then from Definition 12.31 and (|2.6p we deduce 



dt^ 



^m, ^ ={uyB{m)- (2.7) 

at 



From (|2.5p and (|2.7p we obtain the following proposition. 

Proposition 2.5. Let ^ = (^i, . . . ,^fc) he an integrable SOPDE. If ip is an integral 
section of then ip — 4>''^\ where (j)^^'^ is the first prolongation of the map (j) — Toip : 
■^k ^ T^Q Q, and (f> is solution to the system of second order partial differential 
equations 

.^^{t)^{U)W{t).^{t)) l<i<n,\<A,B<k. (2.8) 

Conversely, if (f> : Mf^ Q is any map satisfying i2.8\) then (j)^^"^ is an integral 
section of £, = {^1,. .. ,£,k). ■ 
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From (|2.8p we deduce that if ^ is an integrable SOPDE then {^a)b — fo'^ 
al\A,B = l,...,k. 

The foUowing characterization of the SOPDe's, using the canonical fc-tangent 
structure of T^Q, can be obtained from ((Ti)) and ((T^ . 

Proposition 2.6. A k-vector field = {£^1, . . . on T^Q is a SOPDE if, and only 
= A^, /or an A = l...,fc. 

(Proof) This is a direct consequence of the local expressions (|2.3p . (|2.4p and 
of S''^, and £^a, respectively. ■ 



2.2. Lagrangian formalism: fc-symplectic Lagrangian systems. Consider a 
Lagrangian function L : T^Q ^ M.. We now define the action integral 



(Lo0(i))(t)dH. 



Uo 



where d'^t = dt^ A . . ./\dt^ is a volume form on K^' , : [/q C K.'^ — > Q is a map, with 
compact support, defined on an open set Uq and 0^^^ : J7o C ^ T^Q denotes 
the first prolongation of (p. A map (j) is called an extremal for the above action if 







s=0 



for every flow Ts on Q such that Ts(q) = q for all q in the boundary of (/'(C/q)- Since 
such a flow Tg is generated by a vector field Z G X((5) vanishing on the boundary 
of (/)(C/o), then we conclude that (p is an extremal if and only if 



Uo 



((£z<^L) o0(i)) (t)d'^'t = 0, 



for all Z satisfying the above conditions, where is the complete lift of Z to T^Q. 
d 

Putting Z = Z^—^, taking into account the expression (12. 2p for the complete lift 

Z'^ and integrating by parts we deduce that 4>{t) = (0*(t)) is an extremal of 3 if 
and only if 





"A d 

^ dt^ 

.A=l 




( dL 




dL 


/ 

'Uo 




\dv\ 


0(i)(t)) 





0(l)(t) 



Z'd^t = , 



for all values of Z*. Thus, (j) will be an extremal of 3 if and only if 



A=l 



d 




r dL 




dL 










0(l)(t)) 


" dq" 


0(1) (t) 



(2.9) 



The equations (|2.9p are called the Euler- Lagrange equations for _L. 

We will now give a geometric version of the above equations. 

We introduce a family of 1-forms 6^ on T^Q, I < A < k, using the fc-tangent 
structure, as follows 

et^dLoS"^ 1 < A < fc , 

which are locally given by 

et-^dq^ ■ 

^ dv\ 



If we denote by uof^ — —d9^, in local coordinates we have 



dL 



d^L 



d^L 



8 



M. DE LEON, D. MARTIN DE DIEGO, M. SALGADO, AND S. VILARINO 



Definition 2.7. The Lagrangian L : T^Q — » M is said to be regular if, and only 

( d'^L \ 
if, the matrix is not singular. 

Remark 2.8. Let us observe that L regular if and only if {uj^^ . . . , lu^) is a polysym- 
pletic form and . . . ,lu'1;V — Ker r*), is a fc-symplectic structure (see [35) (See 
appendix for the introduction of some basic concepts on fc-symplectic vector spaces), 
o 

Since {T^Q,uj}^, . . . ,uj^;V) is a fc-symplectic manifold (see Appendix), we can 
define the vector bundle morphism, 

k 

X=(Xi,...,Xfc) ^ r![(Xi,...,Xfc)=trace(zx,wf) = ^ 

A=l 

We now denote by X^(r^iQ) the set of fc-vector fields = (Ci, ■ • ■ , ^l) on T^Q 
which are solutions to the equation 

f7l(ei,...,^i) = d£;i, (2.11) 

where — A(L) — L. The family {T^Q,llj^, El) is called a fc-symplectic La- 
grangian system. If each is locally given by 

d d 

then (Ci, ■ • ■ ,Cl) is a solution to (|2.1ip if, and only if, (^f )* and (Cl)s satisfy the 
system of equations 



dvgdv\ dvgdv\ 
If the Lagrangian is regular, the above equations are equivalent to the equations 

, , ■ v'. + -^-^ (ef )■' = l<z<n, l<A<fc, (2.12) 

and 

i^ty = v\ . (2.13) 

Thus, if L is a regular Lagrangian, we deduce: 

• If = (^1 , . . . , is solution of (PTTTj) . then it is a SOPDE, (see ^A^). 

• There are solutions of (|2.1ip in a neighborhood of each point of T^Q and, 
using a partition of unity, global solutions to (|2.1ip . 

• Since = (Cii ■ ■ • iCl) is a SOPDE, from Proposition 12 . 51 we know that if 
it is integrable then its integral sections are first prolongation 0^^^ : K*^ 
T^Q of maps : R'^ ^ Q, and from we deduce that ((> is solution of 
the Euler-Lagrange equations (|2.9p . 

So, equations (|2.1ip can be considered as a geometric version of the 
Euler-Lagrange field equations. 

• The equation (|2.1ip in the case fc = 1 is t^ujL — dE^, that is the dynamical 
equation of the Lagrangian formalism in Mechanics. 
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Along this paper the family {T^Q,ujf^, El) will be called a k-symplectic La- 
grangian system. 

3. NONHOLONOMIC LAGRANGIAN CLASSICAL FIELD THEORY 

We now bring constraints into the picture. Suppose we have a Lagrangian k- 
symplectic system on T^Q, with a regular Lagrangian L. Let M ^ T^^^Q be a 
submanifold of T^Q of codimension m, representing some external constraints im- 
posed on the system. Although one can consider more general situations, for the 
sake of clarity we will confine ourselves to the case that M projects onto the whole 
of Q, i.e. t(M) = Q and, the restriction t\m ■ M — > Q of t to M is a (not necessarily 
afline) fibre bundle. 

Since M is a submanifold of T^Q, one may always find a covering U of M con- 
sisting in open subsets U of T^Q, with M n C/ ^ 0, such that on each U €U there 
exist m functionally independent smooth functions that locally determine M, 
i.e. 

MnU = {w^ = (wi^, . . . , VkJ e T^Q\ $aK) =0 for 1 < a < m} . 

The assumption that t\j^ is a fibre bundle implies, in particular, that the matrix 
{d^a/dv\){wc^) has maximal rank m at each point G M n J7. 

3.1. The bundle of constraint forms. We now introduce a special subbundle F 
of rank m of the bundle of R'^-valued 1-forms on T^Q defined along the constraint 
submanifold M. The elements rj oi F are K'^-valued 1-forms defined along M which 
are semi-basic, i.e. rj vanishes on the r-vertical vector fields. 

The bundle F is locally generated by m independent E'^-valued 1-forms r]a that 
locally read 

= {vi,---,vt) = ivLdq\---,vt^dq'), (3.1) 
for some smooth functions rj^^ on M C T^Q- The independence of the forms rja 
clearly implies that the m x fcn-matrix whose elements are the functions r|^^, has 
constant maximal rank m (let us observe that m is exactly the codimension of M). 
F is called the bundle of constraints forms. 

Remark 3.1. One interesting case is when F is determined by M through ap- 
plication of a "Chetaev principle" . If the constraint submanifold is giving by the 
vanishing of m functionally independent functions on T^Q, is generated by 
the following R^-'-valued forms: 

77c. = (5^*(rf*a), . . . , ^'=*(rf$a)) = (^ dq\ . . . , ^ dq^) . 

O 

3.2. The constraint distribution. In the sequel, we will show that the constraint 
bundle F gives rise to a distribution 8 along M, called the constraint distribution. 
As above, take F generated by m R'"'-valued 1-forms r/a — {ri\, . . . , 77^) of the form 

(ED). 

Firstly we introduce the following vector bundle morphisms 

Q^r(T,iQ) innnQ) 

X ^ l]^(X) = (7xwi,...,ixu;i), 

where for an arbitrary manifold M we denote by (T^)*M the Whitney sum T*M(B . ^ . 
®T*M of fc-copies of T*M. 
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For each a , {a — 1, . . . , k), let Za G X{T^Q) be the unique local vector field on 
T^Q defined by 

T*(Z„) = and riLiZc) = -Va- (3.2) 

If we write 



d ,„ d 

I 

we deduce from ^JU\) and that 



which determines the {ZaYg uniquely, since L is supposed to be regular. 
One obtains that: 

7 — W'^ 

( d^L \ 

where {W]lg) denotes the inverse matrix of the Hessian matrix . 

The independence of the vector fields Z^ is consequence of the independence of 
the 1-forms rja- Thus the vector fields Za span a m-dimensional distribution §, 
which we will called the constraint distribution. This distribution S will be used in 
Section m 

3.3. The nonholonomic field equations. Summarizing, we are looking for a 
nonholonomic field theory built on the following objects: 

(i) a regular Lagrangian L; 

(ii) a constraint submanifold M ^ T^Q which can be locally represented by 
equations of the form ^ail^ — for a — 1, . . . ,to, where the matrix 
{d^a/dv]^) has maximal rank to; 

(iii) a bundle F of constraint forms and an induced constraint distribution 
S, both defined along M, where F is generated by the to independent 
semibasic R^-'-valued 1-forms p.ip . 

To complete our model for nonholonomic field theory, we now have to specify 
the field equations. We will now introduce the definition of a solution of the non- 
holonomic constrained problem using a generalization of d'Alembert's principle. 

3.3.1. The d'Alembert's principle. Proceeding as in the case of unconstrained field 
theories, we introduce the following definition: 

Definition 3.2. A map (j) : Uq C R'' ^ Q defined on an open set Uq <Z Q with 
compact support, is a solution of the constrained problem under consideration if 
0(1) (C/o) C M and 

f ({Lz^L)o^^^)\{t)dH = Q, 

for each vector field Z on Q that vanish on the boundary of (jjiUo) and .such that 

iZ'V = (3.3) 
for all rj of the bundle F of constraint forms. 

■ d 

Putting Z = Z^—— and taking into account the expression (12. 2p for the complete 



lift Z^, it is easily seen that the condition (13. 3p translates into 
7?^,Z* = 0, l<a<TO, l<A<k, 
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where r]^^ are the coefficients of the constraint forms introduced in (|3.ip . 

One can verify that (j){t) = (t)) is a solution of the constrained problem if 
and only if 





"A d 




f dL 




dL 


/ 






0(1) (t)) 





0(l)(t) 



Z^dJ^t ^ , 



for all values of satisfying (|3.3p . 

Therefore, a solution (p would satisfy the following system of partial differential 
equations: 

fc 

ASr;;^,(0«(t)) (^ = l,...,n) 



dL 
dq' 



0(i)(t) 



y — 



A=l 



/ dL 
[dv\ 



0(1) (t; 



(3.4) 

$a((/.«(t)) = {a = l,...,m). 

As usual, the (a priori) unknown functions play the role "Lagrange multipliers" . 
The equations (|3.4p are called the nonholonomic Lagrangian field equations for the 
constrained problem (compare with A2\ ). 

3.3.2. Geometric description of the nonholonomic Lagrangian field equations. Con- 
sider the following system of equations: 

n{{Xi,...,Xk)-dEL^{n^) , XA^TM,\<A<k (3.5) 

along M. 

One obtains 

Proposition 3.3. Let X = (Xi, . . . ,Xfc) he an integrable k -vector field solution of 
We have 

(i) X = {Xi, . . . , Xk) is a SOPDE. 

(ii) // c/)'"'^-' = {(j)^{t),d(l/ /dt^) is an integral section of X then (j) is solution of 
the nonholonomic Lagrangian field equations \S.4^ . 



(Proof) Let X = [Xi, . . . , Xk) be an integrable /c- vector field solution of (|3.5p . 
Taking into account that rj^ — rj^^dq^ , (see (13. ip ) one obtains that the equation 
p.Sp can be written as follow: 

k 

Y,ix.o^t~dEL^\%V^,dq^ (3.6) 

A=l 

where El = AL — L. 

If each Xa is locally given by 

then [Xi, . . . , Xk) is a solution to (|3.6p if and only if {XaY and {XaYb satisfy the 
system of equations 



d'^L _ d^L 
dq'dv^ dqidv\ 



d^L 
dqidv'^ 



(Xa)' + 



{{Xa)'-v'a) 

d^L 



dL 



dvgdv\ 



d^L 
dv-'gdv\ 



{(XaY v\) = 
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If the Lagrangian is regular, then (Xa)"' — ^^^^ ("'^i' ■ • • ' -^k) is a SOPDE, 
and the above equations are equivalent to the equations 



dL 



(3.7) 



We will now prove (m). 

Let = ((/i'(t),5(/)V9i'^) be an integral section of (Xi, 
a point Wq e M, that is, 



(1) - 



Xk) passing through 
(3.8) 



Substituting 

k 



A=l 



dt^dtB ■ 

in the first group of equations (|3.7p we obtain the equations 



d 




^ dL 


\ dL 











0(l)(t) 



which are the first group of equations in (13.41 
Finally, from 



and since Xa \ jy^- G we have 



{X_-^' 



then (f>QO(/)(i) is a constant function and since 0^^' (0) G M one obtains $c((0(^' (0)) = 
and thus ^a{4''"^H'^)) — 0; that is, the second group of the equations in (|3.4I 
Therefore, we conclude that (f) is solution to the equation 



4. The nonholonomic projector 



The purpose of the present section is to show that for a nonholonomic first-order 
field theory in the sense described above, one can construct, under an appropriate 
additional condition, a projection operator which maps solutions of the equation 
(12. lip for the unconstrained Lagrangian problem into solutions of the nonholonomic 
equations p.5p . 

Given a constrained problem with regular Lagrangian L, constraint manifold 
M C T^Q Slid constraint distribution §, we now impose the following compatibility 
condition: for each Wq G M 

r„^MnS(u;q) = {0}. (4.1) 

If M is locally given by m equations <I'q((?', v'a) — and, if S is locally generated 
by the vector fields Za (see Subsection I3.2p . a straightforward computation shows 
that the compatibility condition is satisfied if and only if 

det(Z„($;3)(u.q)) ^0, 

at each point Wq G M. Indeed, take v G T^j^M Ci §(wq). Then v — v^Zaiwq), for 
some coefficients w". On the other hand, = v{^p) = Za{^(3){wq). Hence, if the 
matrix {Za{^/3)(wq)) is regular, we may conclude that v = and the compatibility 
condition holds. The converse is similar: let us suppose that the compatibility 
condition holds. If the matrix {Za{$f^){wq)) is not regular, then there exist some 
vector V — v^Zaiwq) ^ such that w($/3) ~ and thus v G T^^M n §(wq); 
therefore we conclude that if the compatibility condition holds, then the matrix 
(■^a(*/3)(wq)) is regular. 
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We now have the following result. 

Proposition 4.1. // the compatibility condition (14. ip holds, then at each point 
Wq e M we have the decomposition 

(Proof) The proof immediately follows from (|4.ip and a simple counting of di- 
mensions: 

dimT^^M © §(?i;q) = dimT^u^M + dim §(?i;q) 

— [n + nk — m) + m — n + nk — dim Tu,^ (T^Q) ■ 



We now introduce the following notation: Tj^{TlQ) denotes the restriction of 
T{TlQ) to the submanifold of TlQ, M. 

The direct decomposition of Tj^{T^Q) determines two complementary projection 
operators P and Q: 

P : Tm{TIQ) ^ TM , Q = I-P: Tm{TIQ) ^ S , 

where / is the identity on Tm{T^Q). The projectors P and Q are respectively 
written as follows: 

where (6"'') is the inverse of the matrix {Cap '■ = Za{^i3)). 

The direct sum decomposition of Tjyi{T^Q) determines the following decompo- 
sition of Tl{TlQ) along M: 

Tl{TlQ)^TlM®& 
where for each Wq € M, S^^ is given by 

6^^ =§(wq)® ®§(w;q). 
If § is locally generated by the vector fields Zq, then 

{(Z„,0,...,0), (0,Z„,0,...,0), ...,(0,...,0,Z„), (« = l,...,m)} 

is a local basis of 6. 

The direct sum decomposition of T^{T^Q) along M determines two complemen- 
tary projection operators 7 and Q: 

? : in)M{nQ) ^ nM , Q : [TDMiTlQ) - 6 , 

given by ^(Xi^^, . . . , X^^J = . . . , P(Xfe^J) and Q = 3 - where J is 

the identity on {TI)m{TIQ). Here {TI)m{TIQ) denotes the restriction ofT^{TlQ) 
to the constraint submanifold M. 

Proposition 4.2. Lei = (^]^, . . . ,^^) 6e a solution of the free Lagrangian prob- 
lem, i.e., is solution to the equation ()2.11|) . then M ~ "^{^L ) is a solution 
to the constraint Lagrangian problem. 

(Proof) By definition of T, we know that ^) = (^P{(,1 ■ ■ ■ ' Pi^L ^)) 



with 



e TM. Therefore 



and only if n{ ), . . . , P(ei )) - dE^ £ (ry;^). 



is a solution to (13.51) if 
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We have 



M IM 
fe k 



A=l /l=l A=l 

Thus we conclude that £,l,m = ( P{i\ ),■•■, 5'(fL ) ) is a solution to ((X5 



M 



Remark 4.3. In the particular case fc = 1 we recover the results in [23] o 

5. The nonholonomic momentum equation 

In this section, we derive the nonholonomic momentum equation, the nonholo- 
nomic counterpart to the well-know Noether theorem. More precisely, we prove 
that for every nonholonomic Lagrangian symmetry there exists a certain partial 
differential equation which is satisfied by the solutions of the constrained problem, 
reducing to a conservation law when the constraints are absent. 

Let G be a Lie group and g its Lie algebra. Consider an action G x Q ^ Q. 
The Lie group G acts on T^Q by prolongation of <I>, i.e. 

r,i<i>,K, . . . = ((Tq$g)K), . . . , (Tq<i>,)(z;,J) . 

Definition 5.1. 

(i) We say that the Lagrangian L is invariant under the group action if L is 
invariant under the induced action of G on T^Q. 

(ii) We say that the Lagrangian L is infinitesimally invariant if for any Lie 
algebra element € g we have £,q{L) — 0, where for a vector field X on Q, 
X*^ denotes the complete lift of X of Q to T^Q and is the fundamental 
vector field defined by 

When S,q{L) — 0, then will be called an infinitesimal Lagrangian sym- 
metry. 

Let us now assume that G leaves invariant L, M and F: 

LoTl^g = L, T,}$g(M)cM and {Tl^gr{F)(lF 

for aU 5 e G. 

We consider the vector bundle over Q, defined as follows: denote by 0^(q) 
the linear subspace of g consisting of all ^ G g such that 

CQ(w^q)J^ = for allwq e MnT"^(q). 

We assume that the disjoint union of all g^(q), for all q e Q can be given the 
structure of a vector bundle g^ over Q. 

To any section ^ of g^, one can associate a vector field on Q according to the 
following definition: 

fQ(q):=K(q)]Q(q)- (5.1) 
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Definition 5.2. For each A, the A'^'^ -component of the nonholonomic momentum 
map (J"'')"^ is the map (J"'*)"^ : M (S'ec(g^))* constructed as follows: let ^ he 
any section of , then we define {J^^)~ along M as 

iJ"")f-^^c9t (5.2) 
where is the vector field associated to ^ according to (|5.ip . 

Remark 5.3. In the particular case k = 1, corresponding to the Classical Mechan- 
ics, the above definition coincides with the definition of nonholonomic momentum 
map introduced by Marsden et aHn [8] . o 

Remark 5.4. The map (J"'')^ is the nonholonomic version of the A*'*-component 

^ A 

J(0, . . . , 1^, . . . , 0) ~ (^t^q)) of the momentum map on the polysymplectic man- 
ifolds T^Q defined in [34' when we consider the polysymplectic structure given by 
a;f = -d6lf ,l<A<k/ o 

The relevant role of the nonholonomic momentum map lies in the nonholonomic 
momentum equation. 

Definition 5.5. 

(i) A nonholonomic Lagrangian symmetry is a section ^ of such that 

eH(i) = o. 

(ii) A horizontal nonholonomic symmetry is a constant section of . 

Theorem 5.6. //(/): [/q C R*-' ^ Q is a solution of the nonholonomic field equa- 
tions, then for any nonholonomic Lagrangian symmetry ^ the associated compo- 
nents of the momentum map ( J"'')~ {A — 1, . . . ,k) satisfies the following nonholo- 
nomic momentum equation; 

^ d ^ 

A=l A=l 

along M. 
{Proof) 



^Q^^)-^Qdqr +^A Q^. 



A 

F 



From (|3.4p and taking into account that S,Q{(t>{t)) G 0^ one obtains that the 
above identity is equivalent to 



° ~ dt^ [dv\)^ dqi dv\ ' dt^ dv\ ) [ dtA y ^ dv\ ' ^^-^^ 
where the latter equality is consequence of 



d_ 

ds 



Finally, from ()5.2p and (|5.3p one obtains 

^ dt^ (^"^"'''^f(0(t))) " ^("^"'')^?-(0(t)) 

A=l A=l 
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Corollary 5.7. // ^ is a horizontal nonholonomic symmetry, then the following 
conservation laws holds: 

A=l 

Remark 5.8. If we rewrite this section in the particular case k = 1 we reobtain 
the Section 4.2 in j8j. o 



6. Particular cases 

6.1. Holonomic constraints. 

A distribution Z? on Q of codimension m induces an submanifold M ^ T^Q 
defined as follows: (wiq, . . . , Vkq) is an element of M if VAq G -D(q) for each A{A = 
l,...,k). In coordinates, if the annihilator D'^ is spanned by the 1-forms (pa = 
ipaidq^ (a = 1, . . . , m), then M is the set of solutions to the mk equations = 

If D is integrable, the constraints induced by D are said to be holonomic: in 
this case, takes values in M if and only if (f) takes values in a fixed leaf of the 
foliation induced by £), and we conclude that the constraints can be integrated to 
constraints on Q. 

6.2. Linear constraints induced by distributions on Q. 

The constrained problem: Let Di , . . . , Dk be k distributions on Q and we consider 
the constraint submanifold M = Di © . . . © Dk of T^Q 

If we will assume, for each A {A — 1, ... ,k), that Da is defined by the vanishing 
of niA functionally independent functions on Q, then proceeding as above we 
obtain that the constraint submanifold is given by the vanishing of m = mi + . . . + 
nik independent functions where 

For the bundle F of constraints forms we take the bundle along M, generated 
by the m R'^-valued 1-forms 

i^t^ = J, . . . , = (0, . . . ,rV„,, . . . ,0) . 

Nonholonomic field equations: In this particular C&S6, cl straightforward computa- 
tion shows that the equations p.4p become: 

' k 

V-^f— \-— ^ A"-^ (<^a^),(0(i)) (i = l,...,n), 

^ ^^dt^ t\dv\ 4>'^^Ht)) dq' ^(^Ht) vs^"a7Uv^w; V , , 

^^{(f'^'Ht)) = {aA^l,...,mA,A=l,...,k). 

The constraint submanifold M = D® . ®D. Let D be a distribution on Q. The 
particular case Di — . . . — Dk = D has special interest. As above, if we assume 
that D to be defined by the vanishing of m functionally independent functions 
on Q, then the constraint submanifold M = D® ®D is given by the constraint 
functions 
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We will denote by D'" the distribution on TlQ defined by = {T*ipa) (see 

[m [23] for the case k — 1). Next, we will prove the following two results (see 
Appendix for technical definitions). 

Lemma 6.1. is k-coisotropic in {T^^{T^Q),uj]^{wq), . . . ,uj'l{wq), V{wq)) for 

allw^eM, I.e. {D-)i^cDl^. 

(Proof) In fact, since (D"")*^ is locally generated by semi-basic 1-forms, we deduce 
that 

for all Wq G M, where (D")i^ = {U^,^ e T^Q : u;t{U^^,W^J = 0, V W^^ G 
^w^} denotes the fc-symplectic orthogonal of Dl^^ and Vw^{T^Q) the vertical dis- 
tribution of T^Q at the point w^. ■ 

Proposition 6.2. The following properties are equivalent: 

(i) The compatibility condition holds, that is, TM fl § = {0}. 

(ii) The distribution H = TMDD'" along M is k-symplectic in the k-symplectic 
vector bundle {T{T^Q),luI, . . . ,lu'1,V) 

(Proof) If TM n 5 = {0} then 

TMn§ = TMniD")-^ = 

and 

[TlQ) = T^^M © (D-)^^ Vu;q e M . 
Hence from Lemma |6. II we obtain 

(D-)^^ = (T^^M n {D-)^J © {D'')i^ = H^^ © {D")i^ - © §{wq) . 

Therefore, from a straightforward computation, we obtains that H^,^r\H:^^ = {0} 
or, equivalently, that if^.^ is a fc-symplectic vector subspace of {Tyj^{TlQ),ijj\{wq), 

...,i^i{Wq),V^^). 

Conversely, assume that for each lUq G M, H^, is a fc-symplectic subspace in 
{T^^{TlQ),ujl{wq),...,Lol{wq),V^^), that is H^J^H^^ = {0}. Take 

Z G T^Vi n §(«;q) = T^^M n (i^")i^ c T^^M n (i^'')„^ = il^^ . 

Since ci;^(wq)(Z, y) = for all A (A = 1, ... , fc) and for all Y G -ffiijq, we conclude 
that Z G i?^^. Thus Z G n iJ^^ = {0} and therefore Z = 0. ■ 

Consider now, the restrictions and dnEL to of and dii^L, respectively. 
Since H^^ is fc-symplectic for each Wq G M, there exist a solution on of the 
equation 

k 

J2 ^x.^H - dnEL . (6.1) 

The above equation may be considered as the fc-symplectic version of the char- 
acterization of nonholonomic mechanics in the case of linear constraints given by 

Proposition 6.3. If S,l,m is a solution to the constrained problem, i.e. S,l,m = 
(^i M' ■ ■ ■ ' 'fi 3vt) satisfies (|3.5p . if and only if ^^^^^ is solution of the equation (|6.ip . 
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{Proof) £,L,M is solution of l|3.5p then £,l,m is a SOPDE and jvt ^ T^^^- 
On the other hand, since ^l,m is a SOPDE we have 

along M, and thus Clm ^ (-0")wq- Therefore Clm ^ ^ ^^'^ ^l.m is trivially a 
solution of equation (|6.ip . Conversely, if (Xi, . . . ,Xk) is solution of the equation 
(p?T|) then for each A {A = e C TM and it is evident that 

{Xi, . . . , Xk) is solution of (|3.5p . ■ 

6.3. Linear constraints. 

In the local picture, L is subjected to linear constraints defined by m local 
functions : T^^Q ^ M of the form 

k k 

5=1 B=l 

where //^ be m/c 1-forms on Q, 1 < a < m, 1 < B < k, locally given by /i^ = 

We consider the constraint submanifold 

M^{w^^ K, . . . ,«fej e T^Q : $c.K) = Va} 
of dimension nfc — m. 

We now denote by D the distribution on Q given by 13° = (/if) . 

Proposition 6.4. Let L be a regular Lagrangian and X = {Xi, . . . , Xk) is an 
integrable k-vector field which is a solution of 

k 

Y,ix^u:t-dELe{D^)\ XA\^eTM (6.2) 

A=l 

where (1)^)° =< T*//f >. We have 
i) X = (Xi, . . . , Xk) is a SOPDE. 

a) If (j)^^^ — {(f>^ (t) , d(f>^ / dt^) is an integral section of "K, then <p satisfy the 
equations 



(Proof) Let us observe that -q^ = S^*(d^a) = ^^dq' = ilJ't)idq' = ^Va and 
thus, in this particular case, (_D^)° = {Va)- 



Therefore, the equations (|6.2p are equivalent to equations p.Sp for the case of 
linear constraints. ■ 

6.4. Constraints defined by connections 

Suppose that Q is a fibred manifold over a manifold M, say, p : Q ^ M is a 
surjective submersion. Assume that a connection F in p : Q ^ Af is given such 
that 

TQ^H®Vp, 

where V p = kerTp. We take fibred coordinates (g°, g"), l<a<n — m, l<a< 
m, n = dim Q. The horizontal distribution is locally spanned by the local vector 
fields 
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where stands for the horizontal Hft to Q of a vector field Y on M, and {q'', q^) 
are the Christoffel symbols of F. Thus, we obtain a local basis of vector fields on 

{^.,^«~}- 

Its dual basis of 1-forms is 

{lla = dq\ll^=VZ dq'^ + dq"}. 

We deduce that is locally spanned by the l-forms {i]a}. 
In this situation we have 

T^Q = iJ© .fc. ©i? © Vp ^. ®Vp , 

and for each vector va, 1 < A < k we can write 

o o o o o 

^dq'' ^^(99" "dq""' ^ ^ " ^' dq^" ^ ^ 

We define 

then Wq = (wiq, . . . , Vk^) G M if and only if va^ G H for all A = 1, . . . , /c, which 
means that 

v'X = -v\T'^ foraU A=l,...,fc 

Thus, 

M = {wq e TfciQ : vl = -w^r;^, 1 < A < fc} 
= {wq G Tfcig : - 0, 1 < A < fc} . 

With the 1-forms (/Jq, = F" dq" + dg" on Q we shall consider the 1-forms T*ipa 
on TjjQ. We now consider the equations 

n{{x^,...,Xk)-dEL^{T*^^) 
XaIj^gtm, i<A<fc. 

As in the above example, in this particular situation, these equations are equiv- 
alent to equations p.Sp . 

We write the first equation in (|6.3p as follows 

k m 

Y^ix^ujt-dEL = 51^" 

A=l a=l 

and each Xa as 

From (|6.3p we deduce the three following identities: 



(x^r^i;:!, (x^r-^iS (6.6) 
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If : t/o C K'^ T^Q, V^t) = (V'"(t),?/'"(t),i/'A(t). V'A(t))' is an integral section 
of X = {Xi, . . . ,Xk), then from (j6.4p . (j6.5p and ()6.6p we deduce that 'ip is solution 
to the equations 



VAim) 

d 



k 

A=\ 
k 



dL 
dL 



or, in other words, 

d^p" 

5L 



k 

E 



«S(VXt)) 













dL 




dq^ 











V(t) 



dtp" 



dt^ 

~ — Ar» 



dvj^ 



pa ax 



dL 



^(t) 



dVX dL 



dt dv'X 



These equations are the nonholonomic Euler-Lagrange equations (|3.4p for this 
particular case. 

6.5. The nonholonomic Cosserat rod 

The nonholonomic Cosserat rod is an example of a nonholonomic field theory 
studied in [3T| (see also [H]). It describes the motion of a rod which is constrained 
to roll without sliding on a horizontal surface. 

A Cosserat rod can be visualized as specified by a curve s r{t) in M^, called the 
centerline, to which is attached a frame {di(s), d2(s), d3(s)} called director frame. 
We consider an inextensible Cosserat rod of lenght I. If we denote the centerline 
at time t as s — > r(i, s), inextensibility allows us to assume that the parameter s is 
the arc length. The description of the Cosserat rod can be see in [1T| . 

The nonholonomic second- order model. The model described in [41j fits into the 
multisymplectic framework developed on J^tt, where we have a fiber bundle tt : 
Y ^ X, where usually X plays the role of the space-time and the sections of this 
fiber bundle are the fields of the theory. In this particular case the base space 
X is R X [0, ^] (time and space), with coordinates {t,s) and the total space Y is 
X X X with fibre coordinates {x,y,9). In this model, the fields are the 
coordinates of the centerline (x{t, s), y{t, s)) and the torsion angle 6{t, s) 

Its Lagrangian is given by 



where k = {x"Y + (2/")^i while the constraints are given by 

X + R9y' = and y - ROx' = . 

Here p, a, fi, K and R are real parameters and x = dxjdt., x' dx/ds (analogous for 
y and 9). This model is a mathematical simplification of the real physical problem. 

We now modify this model, by a lowering process to obtain a first-order La- 
grangian: we introduce new variables z = x' and v = y' and obtain the modified 
Lagrangian 



i = f (i' + f) + f - + K{{z'f + (v'f )) + A(z - x') + ^i{v y') 
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where A and fi are Lagrange multiplier associated to the constraint z 



and 



V = y' . This Lagrangian can be thought as a mapping defined on T2Q where 
Q = X S"'^ X = M'^, and if we rewrite this with the notation introduced in the 
section 1 we obtain a /c-symplectic model where the Lagrangian L : Tj^Q ^ M is 
given by 



K 



(6.7) 



subject to constraints 

v\ + Rvlvl = and vl - Rvfvl = . 

In this case the bundle of reaction forces F is generated by the following forms: 
771 = (dq^ + Rv^dq^, 0) and 772 ^ [dq^ - Rv\dq^ ,Q) . 



The nonholonomic fields equations associated to L are given by 



P 



dt^dt^ 



A 



t 5t2 



K 
K 



92 



5*2^*2 



+ 0'(t) 



0'(t) 



9*2 



P- 

R ( A 




0, 



5*1 



9*2 



(6.8) 



where A and /i are Lagrange multipliers associated with the nonholonomic con- 
straints, t = (*^,*^) = (*, s) are the coordinates time and space and the field 
^ : C/q C M2 ^ M'^ are the coordinates of the centerline {(l)^{t),(j>'^{t)) and the tor- 
sion angle (f>^{t). As one can see in the equation ()6.8p the components <j>^, i > A are 
determined by {(j)^ , (fP' , (jy^) . These equations are supplemented by the constraint 
equations (|6.7p . 

Consider the action of R2 x on Q according to the following definition: for 
each (a, 6, 0) e R2 x we consider the map f, e)(9^j ■ • • j 9^) = {q^+a, q^ + b, q^ + 
e,q\...,q'). 

It is easy to see that the following vector field annihilates F along M: 



dq^ dq^ 



This generalized vector field corresponds with the section ^ = (— _R?;|, i?U2, 1) of 



-*r,F 



As S,{L) = the Theorem 15.61 can be applied and the nonholonomic momentum 
equation (see [43]) hence becomes 



R[p 



5201 

Wdt^ ~ 9*2- 



9*2 



R[p 



9*19*1 



90^ \ 901 
9*2" j 9*2" 



9203 



/3 



,2 j,3 



9*19*1 9*29*2 
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This nonholonomic conservation law can also be derived from the nonholonomic 
equations substituting the two first equation in ()6.8p into the three equation. Un- 
fortunately, the knowledge of this nonholonomic conservation law does not help us 
in solving the field equations. 

Remark 6.5. Rewriting the nonholonomic momentum equation in the notation 
used in the description of the second order model, we obtain the nonholonomic 
momentum equation for spacial symmetries given in [53] into the multisymplectic 
setting. o 



7. NON-HOLONOMIC HAMILTONIAN FIELD THEORY 

We now turn to the Hamiltonian description of the nonholonomic system on the 
bundle of fc^-covelocities {T^)*Q of Q. 

The Legendre map FL: T^Q {Tl)*Q is defined (see [T3[2n]) as follows: if 

as s—0 

for each A= 1, . . . ,k and Uq G T^Q. Locally FL is given by 

f)T 

FL{q\v\)^{q\—). 

Assuming the regularity of the Lagrangian, we have that the Lagrangian and 
Hamiltonian formulations are locally equivalent. If we suppose that the Lagrangian 
L is hyperregular, the Legendre transformation is a global diffeomorphism. 

The constraint function on {Tl)*Q becomes = o FL-^ : (T^)*Q M, 
that is, 

dpf 

where the Hamiltonian function H : {T^)*Q — *■ M is defined by H — El ° FL^^. 

dH 

Since locally FL^'^(q\pf ) = {q\ — then 

apf 



H = v\o FL-^ pf -Lo FL-^ 



Thus, from (|3.4p . one obtains 
dp? 



A=l 



where "K^q are the components of the inverse of the matrix (!K^,'^) = (d^H/dpfdp^) 
Note that 



9pf dv, 



c 
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Therefore, the non-holonomic Hamiltonian equations on {T^)*Q can be written 
as follows 



dH 
dH 



^(t) 



^(t) 



-A" 



t 

9*, 



""dpi 



■0(' 



.«gc(*(t))-E|i 



I 



where V : f/ C M'' ^ {TkTQ is locally given by V'(t) = (V''(t), V'j^(t)). 
Next we will give a geometrical description of these equations. 

Let M C {T^)*Q be the image of the constraint submanifold M under the 
Legendre map and let the bundle locally generated by the independent R'^-valued 
1-form 

Va = FL*ria, 1 < a < m . 

Thus, the "Hamilton equations" for the nonholonomic problem can be rewritten 
in intrinsic form as 

k 



E 

A=l 



dHe{ljt), XA\^,eTM 



where rj^ = FL*r]^. 



8. Conclusions 

We have studied various aspects of first-order classical field theories subject 
to nonholonomic constraints in the fc-symplectic framework. The study is very 
similar to the case of particle mechanics (see [24]) and, also, the results are quite 
similar to those obtained in the multisymplectic framework (see (42 J but the k- 
symplectic approach seems simpler in many applications. We have shown that 
in the fc-symplectic approach, the solutions to the equations for the constrained 
problem can be obtained by a projection of the solution to the equations for the 
unconstrained Lagrangian problem. 

We analyze the particular case of a constraint submanifold M which we obtain 
as fc-copies of a distribution on the configuration space Q. In this particular case, 
we construct a distribution H on T^Q along M which is a fc-symplectic subspace 
in {T{T^Q), ujj^, . . . ,uj'l;V) where {ujj^, . . . , uj'l; V) is the fc-symplectic structure ob- 
tained from L. Finally, the nonholonomic momentum map is defined in a similar 
way than in classical mechanics. The applicability of the theory is shown in some 
examples and particular cases. 



Appendix: fc-SYMPLECxic vector spaces 

Let ?7 be a vector space of dimension n(fc-|~ 1), V a subspace of U of codimension 
n and iv^, . . . ,a;'^, k 2-forms on U. For each A {A = 1, . . . , k), ker ui^ denotes the 
subspace associated to uj^ given by 

ker uj^ = {ue U/uj^{u, v) = y v £ U} . 

Definition 8.1. (w^, . . . ,uj''; V) is a k-symplectic structure on U if 

k 

u;^|y^y = 0, f]keiLU^ = 0. 

A=l 
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We say that {U,lo^, . . . V) is a k-symplectic vector space. 

Let ly be a liner subspace of U. The fc-symplectic orthogonal of W is the linear 

subspace of U defined by 



Proposition 8.2. The k-symplectic orthogonal satisfies 

(i) Ac B^B^ c A^. 

(ii) W c (W^)^. 

Remark 8.3. Unlike what happens in the symplectic vector spaces, in our context 
in general dimVI^ + dimW^-"- ^ dimfj. In fact, considering for instance, the real 
space equipped with the 2-symplcctic structure defined by: 

o)^ = A = A V = ker 

where {e^, e^, e"^} is the dual basis of the canonical basis {ei, 62, 63} oi U = R^. 
We consider W = spanjea}, the 2-symplectic orthogonal of W is W-^ = spanjea}. 
In this case, dim W + dim W-^ = 2 ^ dim . o 

We can now introduce the following special types of subspaces of a fc-symplectic 
vector space, generalizing the corresponding notions from symplectic geometry. 

Definition 8.4. Let {U,oj^ , . . . ,oj'';V) be a k-symplectic vector space and W a 

linear subspace of U . 

• W is called isotropic if W C . 

• W is coisotropic ifW^ C W. 

• W is Lagrangian if W = . 

• W is fc-symplectic ifWCi = 0. 

Proposition 8.5. For every vector subspace W of U the following properties are 
equivalent: 

(i) W is an isotropic subspace. 

(ii) v)=0{A=l,...,k) for all u,veW. 

(Proof) Let us suppose that W is isotropic, then W C W-^. Therefore if u, v € 
W then v) = since u G W""-. 

Conversely, u G W, then for each v € W, uj^{u, v) = {A = 1, . . . ,k). Therefore 
u G W-^. Thus we can conclude that W C W-^. ■ 

Proposition 8.6. Let (U,oj^, . . . ,uj^] V) be a k-symplectic vector space and W a 
linear subspace of U . W is a k-symplectic subspace if and only if W with the 
restriction of the k-symplectic structure of U to W is a k-symplectic vector space. 

{Proof) Let us suppose that is a /c-symplectic subspace oi {U,ui^ , . . . ,lv'';V). 
Consider the restriction to W of the 2-forms w^, we will now prove that 
{W, cjIy, . . . , io^, V n W) is a fc-symplectic vector space. 

Given u,v €V nW one obtains 



On the other hand, if u G PI keiu^^ then co^{u, v) — {A = 1, . . . ,k) for all w G 
W, then u G W-^. Therefore, since u G W^fl VF-'- = {0} we deduce nkerwjy = {0}. 

Conversely, ifuGWn W-'- then u G W-'-, that is (j^{u,v) = {A = 1, . . . , k) for 
all V G W. Since u gW we obtain that u G fl ker = {0} an therefore u = 0. ■ 



W-^ = {uG U/ = for all u; G W, A = 1, . . . , fc} . 




uj'^lvxviu^v) ^ {A = 1, . . . ,k) . 
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Definition 8.7. (Awane [S]) A fc-symplectic structure on a manifold M of dimen- 
sion N — n + kn is a family {cu^ , V;l < A < k), where each uj^ is a closed 2-form 
and V is an integrable nk- dimensional distribution on M such that 

(i) w^|vxy = 0, (m) n^^i kerw^ = {0} . 

Then {M,uj^ ,V) is called a fc-symplectic manifold. 

Let us observe that if {M,u!^ , . . . ,uj'' ,V) is a fc-symplectic manifold, then for 
each X G M, we have that (cj^, . . . , uj^, Vx) is a fc-symplectic structure on the vector 
space TxM. 
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